Logarithmic Relaxation in a Colloidal System 
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The slow dynamics for a colloidal suspension of particles interacting with a hard-core repulsion comple- 
mented by a short-ranged attraction is discussed within the frame of mode-coupling theory for ideal glass tran- 
sitions for parameter points near a higher-order glass-transition singularity. The solutions of the equations 
of motion for the density correlation functions are solved for the square-well system in quantitative detail by 
asymptotic expansion using the distance of the three control parameters packing fraction, attraction strength and 
attraction range from their critical values as small parameters. For given wave vectors, distinguished surfaces 
in parameter space are identified where the next-to-leading order contributions for the expansion vanish so that 
the decay functions exhibit a logarithmic decay over large time intervals. For both coherent and tagged particle 
dynamics the leading-order logarithmic decay is accessible in the liquid regime for wave vectors of several times 
the principal peak in the structure factor. The logarithmic decay in the correlation function is manifested in the 
mean-squared displacement as a subdiffusive power law with an exponent varying sensitively with the control 
parameters. Shifting parameters through the distinguished surfaces, the correlation functions and the logarithm 
of the mean-squared displacement considered as functions of the logarithm of the time exhibit a crossover from 
concave to convex behavior, and a similar scenario is obtained when varying the wave vector. 

PACS numbers: 61.20.Lc, 82.70.Dd, 64.70.Pf 



I. INTRODUCTION 



The dynamics in an interacting many particle system is con- 
veniently described by density autocorrelation functions 
for time t and wave vector q. These correlation functions can 
be measured in both experiment and computer simulation 1 1 ]. 
Mode-coupling theory for ideal glass transitions (MCT) dis- 
cusses the transition from a liquid to a glass as a bifurcation 
in the long-time limit of the correlator § q (t) Q]. In the liq- 
uid state, the correlation function decays to zero. If a control 
parameter, say density, exceeds some critical value, the long- 
time limit changes discontinuously from zero to a finite value, 
a glass transition occurs This liquid-glass transition is 
identified with an A%- or fold singularity |4] in the equations 
of motion of MCT. The simplest example for a liquid-glass 
transition is found in the hard-sphere system (HSS), where the 
interaction potential among the particles is zero unless their 
mutual separation becomes smaller than their diameter where 
the potential becomes infinitely repulsive, thus preventing the 
particles from overlapping. The HSS is the system MCT was 
applied to first |3], and it is also the system for which the 
most detailed predictions have been worked out |5, 6]. Close 
to the singularity, the equations of motion can be expanded 
in asymptotic series. This yields a two-step decay with two 
related power laws for the short-time and the long-time de- 
cay at a liquid-glass transition |2]. The HSS can be realized 
in colloidal suspensions |7]. Experiments in these systems 
lead to the conclusion that MCT is able to describe the main 
aspects of the glass transition qualitatively and some aspects 
even quantitatively 0] • 

MCT can also exhibit other singularities than the fold |2]. 
These higher-order singularities were predicted recently to oc- 
cur for colloidal systems with short-ranged attraction where 
A3- and A4-singularities were found that are also called cusp 
and swallowtail I9I [Tol [Till . In these systems, the hard-core 
repulsion is supplemented by a short ranged attraction, e.g. 



in the square-well system (SWS). A cusp singularity is the 
endpoint of a line of glass-glass transitions that arises if two 
different mechanisms of arrest are of the same importance. In 
the SWS the first mechanism is the hard core repulsion that 
causes a transition as in the HSS via the well known cage ef- 
fect. The second mechanism leading to arrest is bond forma- 
tion introduced by the attractive part of the potential. This 
latter transition was proposed as relevant for the transition to 
a gel 1 10]. If the difference in the two mechanisms is less pro- 
nounced, the glass-glass transitions vanish and give rise to an 
A/t-singularity. In the SWS this happens as the range of the at- 
traction is increased [11]. The range of attraction considered 
here is of order less than 20% of the particle diameter and the 
strength is about several k^T. This is accessible in colloid- 
polymer mixtures with nonadsorbing polymer which is well 
under control experimentally 1 12]. Higher-order singularities 
have also been identified for a number of short-ranged po- 
tentials with shapes differring from the SWS yielding certain 
quantitative trends but no qualitative changes 1 13]. 

In addition to the success of MCT for the description of 
the HSS, two findings support the use of this theory for the 
description of colloids with attraction. First, a reentry phe- 
nomenon was predicted by the theory where a glass state is 
melted upon increasing the attraction @, [Till . This was sub- 
sequently found in several experiments |l4 fl^l and com- 
puter simulation studies fl5l flrj fl7[ Il8ll. Second, there 
are indications of logarithmic decay fl9ll and related anoma- 
lous decays 1171 .12011 that are consistent with scenarios found 
within MCT l^ fTlll . To investigate the dynamics in such sys- 
tems, apart from computer simulation dynamic light scatter- 
ing has already been used to determine correlation functions 
I14fl5ll20ll . Direct imaging techniques are available to deter- 
mine also the mean-squared displacement (MSD) with high 
precision |2l|, |22j, 12311 . The purpose of the present paper is 
the application of the general theory for higher-order singu- 
larities 1 24] to the SWS and the derivation of testable quanti- 
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tative predictions for the correlation functions and the MSD. 
Certain scenarios have been discussed before for schematic 
models l24ll . Some of these scenarios shall be identified also 
in the microscopic model in the following. 

The paper is organized as follows. In Sec.[n]the equations 
of motion and the asymptotic solution for the logarithmic de- 
cay are summarized and the subdiffusive power law for the 
MSD is derived. The theory is applied to the A4-singularity in 
the SWS for the correlation functions in Sec. IIIII and for the 
MSD in Sec. lIVI Changes in the scenarios when moving from 
the A4-singularity to A3 -singularities are discussed in SecM 
and Sec. lVII contains a comparison to results obtained for the 
hard-core Yukawa system (HCY). Sectio rfVIll presents a con- 
clusion. 



II. ASYMPTOTIC SOLUTIONS 

We shall consider a system of N particles with diameter d 
in a volume V interacting with a spherical potential. When at 
time t the j'th particle is located at fj{t) the density variables 
are defined as p q (t) = £ ; exp[/^r ; (f)]. 



A. Equations of Motion 

The equations of motion for the normalized density corre- 
lators ty q (t) = (p|(*)P?) /(\Pq\ 2 ) within MCT, when Brownian 
dynamics for the motion in colloids is assumed, are given by 

iSSillli, 



x q d t ty g (t)+ty q (t)+ J m q (t-t')d t ,i! q (t')dt' ■■ 



0. 



(la) 



The initial conditions are 0^(0) = 1. The microscopic time 
scale reads i q = S q / (Doq 2 ). It is given by the short-time diffu- 
sion coefficient, Do, characterizing the Brownian motion and 
the static structure factor S q — (|p^| 2 )- The mode-coupling 
approximation results in expressing the kernels m q (t) in terms 
of the correlators ty q (t) |01, 



(lb) 



As a consequence of the factorization into pair modes for the 
structural relaxation in simple liquids, f q is a bilinear func- 
tional of the density correlators [3], 



- 1 f d 3 k , ~ 



(lc) 



and the vertex is determined completely by the static structure 
of the liquid system 0H, 



v aM = s <i s k s \a-k\P q- kc k+q-(q- k ) c \ q -k\ 14 



(Id) 



The number density is given by p = N/V and c q denotes the 
direct correlation function which is related to the static struc- 
ture factor S q in the Ornstein-Zernike relation, S q — 1/[1 — 
p Cq] , both depend on external control parameters like density 



or temperature |29]. For the SWS with hard-core diameter 
d, depth of the potential uq, and range of the potential A, we 
get three dimensionless control parameters, the packing frac- 
tion cp = d i pn/6, the attraction strength T — UQ/{k^T) and 
the relative well width 8 = A/d. These can be combined to a 
control-parameter vector V = (cp,r,8). 

It is the long-time limit of the correlation function, 
lim f _>oo(|) 9 (f) = f q , that determines whether a system is in the 
liquid regime, where f q = 0, or in an arrested state, where 
< f q ^ 1. In the latter case, the values f q characterize the 
arrested glassy state and the f q are called glass-form factors 
or Debye- Waller factors. In the long-time limit, the equation 
of motion Eq. Q reduces to an equation involving only the 
mode-coupling functional and the glass-form factors 0], 



(2) 



Frequently studied is the dynamics of a single or tagged 
particle with the single particle density (?) = exp[iqr s (t)]. 
For the correlation function of a tagged particle, § q (t) = 
(pj£(f)p|), similar equations as Eqs. {lj have been derived 



^3^(0+^(0 + J q >n q (t - t')M q {t')dt' = 0, (3a) 

with v q — 1/ (D s Q q 2 ). The short-time diffusion coefficient for 
a single particle, D , again specifies the Brownian dynamics. 
The mode-coupling functional for the tagged particle motion, 



d 3 k 



(27t) 3 q- 



S k ^cf(qk) 2 f k f , , (3b) 



is also determined by the static structure of the liquid system 
where c s q is the single -particle direct correlation function 12911 . 

The dynamics of the tagged particle is coupled to the co- 
herent density correlator §q(t) and for that reason (j)^(f) also 
displays the bifurcation dynamics that is driven by (j) 9 (f ). The 
equation for the long-time limits of the tagged particle corre- 
lations function, (^(f — * °°) = f s reads 



/|/(W|) = ^LA/ S ]- 



(4) 



In the following, the tagged particle will be assumed as of the 
same sort as the host fluid. If the host particles are in the liquid 
state, f q = 0, a tagged particle cannot be arrested, and in that 
case Eq. @ implies f q = 0. 

The MSD is defined by 5r 2 (r) = {\r s (t) -r s (0)\ 2 ) and de- 
scribes the average distance a particle has traveled within 
some time / l29Tl . It is obtained, e.g., as small wave-number 
limit of the tagged-particle correlator in Eq. yj, § q (t) = 
\-q 2 hr 2 (t)/6+0{q A )[Mk, 

5r 2 (t)+D s [ mPXt-t'^it^dt' = 6D s t, (5a) 



o 

m m (t) = lim ? ^ «4(f) = T MSD [W)A S (0]- The mode- 
coupling functional for the MSD reads 



TMSD[f,f] 



(5b) 
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A characteristic localization length r s is defined by the sec- 
ond moment for the relaxation of the distribution of § s a (t) 
01, which can be identified with the functional in Eq. J5bi 
r 2 = 1 / TMSD[fJ s ]- It is the long-time limit of the MSD. Its 
value at the critical point, r£, characterizes the arrested struc- 
ture. The value 6r c s 2 represents the plateau for the dynamics 
of8r 2 (f). 

Equations Q to (O are solved numerically using algo- 
rithms introduced in 1 30, 31]. Details of the implementation 
are found in |5|, [Till . We use d as unit of length, d = 1 , and 
choose the unit of time so that 1/Dq = 1/Dq= 160. The struc- 
ture factors for the SWS and the HCY are calculated in mean- 
spherical approximation 1 1 ll 13211 . The wave numbers shall be 
discretized to a grid of M points with a spacing Aq = 0.4/ d. 
The cutoff in the calculations is ranging from M = 300 for 
8 > 0.04 up to M = 750 for 8 = 0.02. 



B. Logarithmic Decay Laws 

The asymptotic solution at higher-order singularities shall 
be quoted from Ref. 1 24] where also further details can be 
found. The asymptotic expansion is performed in small devi- 
ations of the correlation function from the critical long-time 
limit /' involving the coefficients 



qk 



which can be split into values at the singularity, A^ c _ k , 



(6) 



(n) 

and remainders, A y J,„ k (V) = A 



An)c 



r(») 



qkykn' 

./M qkyk^^qkykS^)- TheJa " 

cobian matrix of Eq. (|2) is singular at the critical points and 
assumes the form [8^ — A^ c ] . The non-negative left and 

right eigenvectors of matrix A^ L shall be denoted by a* and 
a q and can be fixed uniquely by requiring Y. q a* a q = 1 and 

Y,q a q a q = 1- The reduced resolvent R q k of A^ c maps vectors 
orthogonal to a* to vectors orthogonal to a q . The leading- 
order result for Eq. Q is then given by 



$ q {t)=f c q +h q [-B\n{t/%)] , B = y [— 6ei (V)/tc 2 ] , (7) 

with the critical amplitudes h q — (1 —f q )a q and the separation 

parameter £i(V) = a£Ag (V), which is restricted to negative 
values, Ei < 0. If e indicates the distance of the control param- 
eters V from the critical point, the leading result is of order ^/i 
and correct up to O(e). The next-to-leading-order approxima- 
tion is 

<M0 = (£+/*) + h q [{-B+B x )\n{t/^) 

+ {B 2 +K q B 2 )\n 2 (t/z) 

+B 3 ln 3 (fA)+B 4 ln 4 (f/T)] , (8) 

what includes the terms of order e and neglects terms of order 
e 3 / 2 . It involves corrections to the plateau values, 



f q = {l-f q )R qk 



A[ 0) (V)- £l (V) fl 



(9) 



correction amplitudes, 



K q — R q k 



and the prefactors, 



4 (2)c 2 
A kkik 2 a ki a k 2 a k 



/a q , 



(10) 



B\ = (0.44425 C- 0.065381 W )ei(V)- 0.22213 e 2 (V) , 

(11a) 



B 2 = (0.91189 i; + 0.068713/j 3 )ei(V)-0.15198e 2 (V) , 

(lib) 



Bo, = -0.1 3504^3 6i(V) , B 4 = -0.046197 W Ei (V) . 

(11c) 

The numbers characterizing the higher-order singularities are 



^ = Y j a;[a 2 K q + a q /2] 



(12) 



and 



.(3)c 
qk\k 2 k 3 



ak,ak 2 a k 3 +lA q 2 k^k 2 a ki a k 2 Kk 



(13) 

For the leading correction also an additional separation pa- 
rameter is introduced, 



£2(V)=I«S ) (VK + 2E,(V)j;fl^ 

q q 



2 L< A % 2 akjk 2 /{1 -f h ) ~ a q f q /{\ -ft 



(14) 



The time scale x is determined by matching asymptotic ap- 
proximation and numerical solution of tyq(t) at the plateau f q 
or the rescaled plateau fjj + f q for Eq. Q and Eq. ©, respec- 
tively. 



C. Coupled Variables 

Inserting the asymptotic expansion of Eq. into the long- 
time limit of Eq. d3al . the approximation for the tagged parti- 
cle dynamics up to order e is obtained 13311. 

+ (B 2 +K°B 2 )ln 2 (t/%) 
+B 3 ln 3 (f/x)+B 4 ln 4 (f/x)] , (15) 

with the critical amplitudes h s q = (1 — f q c )a q , the correc- 
tion amplitudes K q and the plateau corrections fi. The latter 
are derived from the functional Obt and the related coherent 
quantities by 



(16a) 



4 



it k 

+ Y} A q,kp a k aS p + A qk P aka P + A qk,p Clka p 1 \ ' 



(16b) 



£(S # -A^)a£/* = -ei(V) fl f +£a^/*+A^(V). 

(16c) 

The derivatives with respect to the coherent and tagged parti- 
cle glass-form factors are denoted before and after the comma, 
respectively. The coefficients are 

A s m 1 Lq r) g^lM^S 



x(l-/ fc e 1 )-(l-/^)(l-/A e )-(l-/; n c ) = 

= 4 se -(-.A* (V) 



Similar arguments as above yield the asymptotic expansion 
for the MSD up to order e , 

l -hr 2 {t) = r c2 - r 2 s -hASD[{-B+Bi)\n{t/x) 

+ (B 2 +K MSD B 2 )\n 2 (t/z) 
+B 3 ln 3 (f/T)+B 4 ln 4 (f/T)] , (18) 

with parameters 

huso = r c s 4 {T^ SD [h k ,f^} + TfaMtfp]} , (19a) 



K M SD =r C s {fMSD[h,K] + ^MSdIM^/T] 



+ ?MSD[/fc, W]}/%SD ~ AmSd/1* 



.c2 



(19b) 



to account for the fact that the MSD is conveniently shown 
in a double-logarithmic representation which is more sensi- 
tive to the detection of power laws. The asymptotic approx- 
imation (1181 for the MSD can be written as z = «o + d\y + 
a 2 y 2 +a , iy i +a4y 4 . Here, z = dr 2 (t)/6 and y = ln(f/x) The 
constant term represents the square of the corrected localiza- 
tion length, flo = r c s 2 — r 2 , the coefficients a\ = /imsdCB — B\), 
a 2 = -h MSD (B 2 +K MSD B 2 ) as well as a 3 = -/imsd#3 and 
«4 — IimsdB4 are the separation dependent prefactors for the 
leading and next-to-leading order terms. This yields the ex- 
pansion 

lnz = Inrf - r 2 Jr c2 +x' 'y + b 2 y 2 +a 3 /rfy 3 +a 4 /r c2 y 4 + O (e 3 / 2 ) , 

(20a) 

with 



x =a\jr c2 , b 2 



2r L s 2 a 2 — a 2 
2^ 



^\ In leading order, one gets a power law for the MSD, 



Sr*(0/6 = rf {t/xf, 



with an exponent 



x = h M sT>B/r c s 



c2 



(20b) 



(21a) 



(21b) 



Exponent x varies with the square-root in the separation pa- 
rameter £1, cf. Eq. @. Including the corrections of order e 
rescales the exponent to 



x'=h mD {B-B 1 )/r c2 . 
and the next-to-leading order result reads 

8r 2 (f)/6 = (t/%f{r? -? 2 + b 2 rf ]n(t/xf 

ln(f/T) 3 +a 4 ln(f/T) 4 }. 



(22a) 



(22b) 



=^ 4 {^msd[%A./; c ] + fmM^ s J S p\ 

+ ?uMJpW) - ^sd[/(,/; c ](V c ')}Amsd 
-et(V)4 SD /rf . 

(19c) 

For the generic liquid-glass transition, the asymptotic expan- 
sion was carried out with a different convention as in Eq. 
however, the quantities h q , K q , hi,Ki, Ii^sd, ^msd> an d C 
are the same in both descriptions |5l|6J|. The plateau correc- 
tions, f q , fq, and r 2 are different for liquid-glass transitions 
and higher-order singularities. The expansions in Eqs. 0, 
dl 51 , and dl 81 share the coefficients B, B\, B 2 , B$, and B4. 
They differ in the plateau and its correction, the critical am- 
plitude and the correction amplitude. 

D. Subdiffusive Power Law in the MSD 

The logarithmic decay laws shall be phrased for the MSD in 
a slightly different form than in Eq. Jl 81 . This is done in order 



III. CORRELATION FUNCTIONS NEAR AN 
A4-SINGULARITY 

Before we can apply the asymptotic expansion of Eq. (jSJl, 
we need to specify the values for /J3 and £ appearing in the 
prefactors of Eqs. dl It . The A4-singularity is characterized 
by /J3 = 0. This condition has been used to locate the A4- 
singularity at V = V* for the SWS by: 

cp* = 0.52768 , T* = 4.4759 , 8* = 0.0438 1 . (23) 

The vanishing parameter /J3 implies a considerable simplifica- 
tion in the preceding formulas since B3 = B4 = 1 24| . The 
deviations of the control parameter values specifying the A4- 
singularity from the ones reported in Ref. full originate from 
refined numerical procedures used here and they do not ex- 
ceed 6%. The characteristic parameter was < 5 • 10~ 4 at 
the control-parameter values specified above. The parameter 
^ varies regularly and is ^ = 0. 122 at the A4-singularity. This 
is smaller than the value in the HSS, ^hss = 0.269 |5]. 
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FIG. 1: Wave-vector dependent amplitudes characterizing the A4- 
singularity, Eq. 1231 . for coherent and tagged particle correlators of 
the square-well system (SWS). In the upper and middle panel the 
critical glass-form factors /*, Eq. J2j, and the amplitudes h* are 
shown as full lines, respectively. The dashed lines represent the val- 
ues for f q *, Eq. J4j, and h", Eq. <16at . For the hard-sphere system 
(HSS), f q c and h s q are shown dotted. The lower panel shows the 
correction amplitudes K*, Eq. <10> . and K", Eq. J16b> . as full and 
dashed lines, respectively. A square at q = 24.2 indicates the cor- 
rections for the path calculated for Fig. [2] and the correlators shown 
in Fig. [3] The correction amplitudes K q (— ■ — ) and K s q (■■■) for the 
HSS are shown for comparison. The unit of length here and in the 
following figures is the hard-core diameter d of the particles. 



The second prerequisite for the asymptotic description ac- 
cording to Eq. (jSJi are the wave-vector dependent amplitudes 
ft, h* and K*. These are shown for the A4-singularity in Fig. 
[^together with the related values for the tagged-particle cor- 
relator, Eq. (I15> . The quantities for the tagged particle motion 
are close to the ones for the coherent correlator § q (t) except 
for values of q smaller than, say, q =10. This difference was 
observed already for the HSS |6|. Since we will not be con- 
cerned with small q in the following, we restrict the discussion 
to the coherent dynamics and imply that the same is applica- 
ble also to the incoherent part with only minor changes. In 
comparison to the HSS the /*, h*, f q * and h s * are extended 
over a broader g-range. The maximum in hi is shifted from 
q ps 13 to q w 25 reflecting the smaller localization length in 
the SWS as noticed before, cf . [11]. We see in the lower panel 
of Fig.^that the distributions of the correction amplitudes K q 
and K q share that trend of becoming broader from the HSS 



to the Azt-singularity of the SWS. The zero in K q moves from 
around q ps 14 in the HSS to q pa 32 in the SWS. In addition, 
the amplitudes are shifted to lower values for small q. 
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FIG. 2: Curves of vanishing quadratic correction in Eq. {5} at the A4- 
singularity of the SWS, B 2 (q) = (dash-dotted), for q = 7.0, 20.2, 
24.2, 27.0, and for K q = as labeled. The full line shows a part of the 
glass-transition diagram for constant 8 = 8*. The lines of vanishing 
separation parameters Ej (V) and £2(V) are shown by a broken and a 
dotted line, respectively. For the wave vector q = 24.2, a path on the 
curve S2(24.2) = is marked (+) and labeled by n, for which the 
correlators are shown in Fig. [3] State n = 2 is analyzed also in Fig. 
|4| For the points (A) labeled a, b, and c the decay is shown in Figs. [5] 
and[l2l 

Having specified the characteristic parameters for the A4- 
singularity, the solution at any point in the control-parameter 
space can be compared to the asymptotic approximation as 
the control parameters are translated into separation parame- 
ters £1 and 62. As done for the schematic models in Ref. 12411 . 
we start by determining the surfaces where the quadratic cor- 
rections in Eq. ^ are zero, #2(9) = B2 + K q B 2 = 0. On these 
surfaces in the control-parameter space the logarithmic decay 
is expected to show up as straight line around the plateau /*, 
as the cubic and quartic terms in Eq. ^ vanish because of 
Z?3 = B4 = /i3 = at the A4-singularity, cf. Eq. il Id . We get 
a different surface for each wave vector q and show typical ex- 
amples in Fig.|2]for a cut through the glass-transition diagram 
for 8 = 8*. For q = 7.0 one gets K q = — 1 .8 1 . The solution of 
B 2 {\) = 1.81B(V) 2 yields the chain line labeled B 2 {7.0) = 
in Fig. [2] and is lying in the arrested region close to the line 
of liquid-glass transitions. Since the K q depend smoothly on 
q, the evolution of the curve where B 2 (q) = 0, can be under- 
stood by inspecting the parameters B and B 2 . The square B 2 is 
always positive and proportional to 61 (V), cf. Eq. Q, there- 
fore K q B 2 is proportional to K q (V) | and shares the sign of 
K q . Inserting /J3 = and £ = 0.1216 into Eq. ( II Vol yields 
B 2 (\) =0.111ei(V)-0.152e 2 (V), which has to be positive 
to comply with B 2 (q) = 0. The second separation parameter 
is negative, z 2 (V) < 0, below the dotted curve for z 2 = in 
Fig. ED In addition, the value |e 2 (V)| on the line B 2 (7.0) = 
is larger than |£i(V)| which we can also infer from the fact 
that the line £1 = is closer than the line e 2 = 0. We now 
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chose a point on the line Z?2(7.0) = 0, keep the first separation 
parameter fixed, say Ei = e'j, and move to higher values for 
K„, e.g., for q = 20.2 where K q = —0.966. B 2 stays the same 
and the term K q B 2 increases. To ensure that Z?2(20.2) = 0, the 
value B2CV) has to decrease. We can achieve that by mov- 
ing closer to the line £2 = 0. For fixed e\ this implies a shift 
to lower cp and higher F. Consequently the entire line, where 
Biiq) — 0, is rotating clockwise around the A4-singularity as 
K q increases. This is seen for the chain line Z?2(20.2) = in 
Fig. I2 Since K q is monotonic increasing with q and Kq has 
the same trend when neglecting the small oscillations, Fig.^ 
the line B s 2 (q) = B 2 (V) + B(\) 2 K s q = also rotates clockwise 
with increasing wave-vector q. 

The variation of the lines B2(q) — described above de- 
pends only on the angle at which the lines £1 = and £2=0 
intersect at the Azt-singularity. This intersection is in a sense 
generic that it is shared by the close-by A3 -singularities of the 
SWS. It applies also to the A/t-singularities of the other poten- 
tials which are similar to the square well. This is so because 
the functionals determining the separation parameters depend 
on quantities like the structure factors and the glass-form fac- 
tors which are similar for different potentials 1 1 3 ] . For a given 
wave vector q, the line B2 (q) — may or may not lie in the liq- 
uid regime depending on K q . For the SWS at 8 = 8* we get a 
range of — 1 < K q < 0.4 corresponding to 20 < q < 35, where 
a line B2 (q) = is found in the liquid regime. We illustrate 
this by adding lines for q = 24.2, q = 21.0 and for K q = to 
Fig-El The vanishing K q is corresponding to q m 32.3 yielding 
a line B2(q) — still in the liquid. For q > 35, the latter line 
rotates further around the A4-singularity and into the arrested 
regime beyond the almost horizontal line of liquid-glass tran- 
sitions. 

We select a wave vector q = 24.2 with K q = —0.596 as in- 
dicated in Fig. [2 by a square and choose a path on the line 
#2(24.2) = marked by the plus symbols in Fig. |2] For n = 
1,2,3, the control parameters are (r,<p) = (3.312,0.5125), 
(4.271,0.5250), and (4.453,0.5274), respectively. The solu- 
tions are shown in Fig. [3] together with the leading approxi- 
mation, Eq. Q, (dotted) and the next-to-leading approxima- 
tion, Eq. (|8j, (dashed). The time scales x are matched at the 
plateau f* for the leading approximation and at the renormal- 

ized plateaus f q t +f q for the first correction. We recognize that 
for n = 3, Eq. i|8} accounts for more than ten decades in time 
with a relative accuracy better than 5%. The leading approxi- 
mation is acceptable on that level for nine decades. For n = 1 
two and more than one decade are covered, respectively. Five 
and three orders of magnitude in time are achieved for n = 2. 
For n — 1, 2, 3, the leading approximation describes at least 
30% of the complete decay and when including the correction, 
65% are covered on the chosen accuracy level. The distance 
in the control parameter Y from the value at the A4-singularity 
is 25% for n = 1 and 4% for n = 2, so no fine-tuning was nec- 
essary to obtain such large windows for the logarithmic decay. 
The curve n = 1 requires about five decades for the complete 
decay which is well in the reach of today's computer simula- 
tion techniques 113411 . 

It was possible to describe part of the critical decay at an 
A3-singularity in a one-component model by the expansion in 




FIG. 3: Logarithmic decay at the A4-singularity in the SWS for 
q = 24.2 on the path indicated in Fig.|2| The correlation functions are 
shown as full lines for the states n = 1, 2, 3 (see text) and at V = V* 
labeled <K. The horizontal line indicates the critical plateau value /* 
for q — 24.2, short lines the renormalized plateaus /* + f q . Broken 
lines show the approximation of Eq. — (B — By) ln(f/x), dotted 
lines the approximation by Eq. 0. Filled and open symbols, respec- 
tively, mark the points where the approximations deviate by 5% from 
the solution. The cross indicates the time when the solution for n = 3 
and the critical correlator 0* differ by 5%. The unit of time here and 
in the following figures is given by a short-time diffusion coefficient 
ofZ> = 1/160. 



polynomials in lnf at a point away from the singularity 1 24]. 
We therefore compare the critical decay (])*(f) with the de- 
cay for n = 3 and indicate the point at f 5000 where both 
differ by 5% in Fig. |3 With only the leading correction at 
our disposal, a 2%-criterion was not fulfilled as for the one- 
component model, where also the next-to-leading correction 
could be used |24]. The dashed line for n = 3 does not come 
closer to the critical decay than 4%. Allowing for 5%, the in- 
terval from f « 20 to f w 4000 could be described. However, at 
the A4-singularity the approximation in Eq. (|8j always yields 
a straight lnf-decay as approximation on the chosen path with 
B2{q) = 0. This disagrees qualitatively with the observed crit- 
ical decay. 

To identify correctly some decay that is linear in the § q (t ) 
versus logf diagram with the logarithmic decay predicted by 
the asymptotic laws, Eq. l|8), we check if a different correlator 
with a different correction amplitude K q is not linear in lnf 
at the same point in the control-parameter space. For a two- 
component model a characteristic alternation of concave, lin- 
ear and convex decay in lnf was found 1 241 13311 . Not both cor- 
relators could be linear in lnf at the same time. For the SWS 
this check is performed at the point n = 2 from Fig.|2|by vari- 
ation of the wave vector q. For the wave vectors q = 4.2 and 
32.2 the correction amplitudes are K q = —1 .400 and —0.0413, 
respectively. Therefore B 2 (4.2) < andfl 2 (32.2) > 0. We ex- 
pect 9 (f) to be concave or convex, accordingly, as is demon- 
strated by the inset of Fig.|4] The rescaled correlators § q (t) 
displayed in the full panel allow for a more detailed analy- 
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FIG. 4: Logarithmic decay at the A4-singularity for varying wave 
vector q. The inset shows the correlation functions tyq(t) at state 
n = 2 from Fig. [2] for wave vectors q = 4.2, 24.2, and 32.2 from 
top to bottom and the short horizontal lines show the corresponding 
critical plateau values ft. The full panel shows the same correlation 
functions rescaled according to <j> ? (f) = (tyq(t) — fq — fq)/h* q as full 
lines and labeled by the respective wave vectors. Dashed lines show 
the asymptotic laws, Eq. {8j. The deviations of the approximations 
from the solutions of 5% are marked by the open symbols. Filled 
symbols for q = 4.2 (T) and q = 32.2 (■) show the 5% deviation 
from the additional approximation of neglecting quadratic terms in 
Eq. {8} (see text). 



sis. We see that the solutions as well as the approximations 
clearly exhibit increased curvature for larger q. Since the 
coefficient linear in lnf is not depending on q, cf. Eq. {8j, 
the middle dashed line represents the leading correction to all 
three correlators when the quadratic terms are neglected. For 
q = 24.2 we observe good agreement over almost 5 decades 
as before, cf. Fig. [3] For q = 4.2 and 32.2, however, the ad- 
ditional approximation reduces the range of applicability to 
less than one decade as marked by the filled symbols. Includ- 
ing the quadratic terms from the approximation (JSJl extends 
this range by half a decade to later times and to earlier times 
by one and almost two decades for q = 4.2 and 32.2, respec- 
tively. The time window defined by a 5% deviation from the 
approximation l|8} is larger by one and two orders of magni- 
tude for q — 24.2 than for q — 32.2 and q — 4.2, respectively, 
what indicates that ^-dependent higher-order corrections sig- 
nificantly influence the range of applicability for the leading 
correction ((SJ. 

The time scale x in Fig.0]was matched for q = 24.2, so the 
violation of scale universality inherent to an approximation 
like in Eq. © leads to different times x(q), where the correla- 
tors for different q cross their respective renormalized plateau 
ft +f q 1 24]. The representation with the rescaled § q (t) is par- 
ticularly sensitive to these deviations since the point where the 
plateau is crossed is required to be zero, <p ? (Y/x) = 0. In Fig. 
|4]we see that the line crossing the zero is slightly broader than 
a single curve. The deviations in x(q) are small enough to 
not exceed the numerical grid for the time axis which around 
X = 2988 is given by At = 172. So we interpolate to get for 



q = 4.2, 24.2, and 32.2, i{q) = 2899, 2988, and 3017, respec- 
tively. These differences do not introduce larger errors in the 
analysis carried out above. 
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FIG. 5: Logarithmic decay at the A4-singularity for the three states 
marked by triangles in Fig. [5] The inset shows the correlation func- 
tions §q(t) for q = 24.2. The plateau value /* is indicated by the 

short horizontal line. The full panel shows = (§q(t) — fq — 
fq)/h q divided by the respective values for (B — B\) at the three 
states specified. The dashed curves show the result from Eq. {8j. 
Filled squares and circles mark the points where curve a and c de- 
viate by 5% from — ln(f/x), respectively. The deviation for curve b 
(A) for short times is at t/x ~ 10~ 4 and not included in the figure. 

In order to change from convex to concave behavior we 
can also change the control parameters. For states above the 
line Z?2(24.2) = 0, we expect concave behavior, B2(q) < 0, for 
states below, convex decay, Z?2(24.2) > 0. For a demonstra- 
tion of this result, the rescaled correlators § g (t/%) at the states 
labeled a, b, c in Fig.|2]are divided by the prefactor (B — B\) 
of Eq. (|8j. This way the part of the decay that is linear in 
In? shows up as straight line with slope — lnlO in Fig. [5] The 
approximations (JSJl are shown as dashed lines for each state 
representing -ln(f/x) + [B 2 (24.2)/(B - Bi )} ln 2 (f/x). For 
state b the approximation is identical to — ln(f/x) and the so- 
lution follows that line over 5 decades before 5% deviation is 
reached. The states a and c are chosen to have the same value 
for B B { w 0.015 and B 2 (24.2) = ^0.0020, respectively. 
The solutions at state a and c follow the —lnf -law closely 
within a 5% margin for two decades or one decade, respec- 
tively, which is significantly less than found for state b. We 
can infer from Fig. [2] that at state a the quadratic corrections 
would vanish again if we went from q = 24.2 to the higher 
wave vector q = 27.0. A scenario similar to the one shown in 
Fig-Hcan be found. 

The procedure outlined in Figs. EI and|5]can be sum- 
marized as follows. From the higher-order singularities there 
emanate surfaces in the control-parameter space for a specific 
wave vector q where the quadratic term in Eq. ^ is zero, cf. 
Fig. |2j and the decay is linear in Inf. Moving closer to the 
singularity on that surface, the window in time where the log- 
arithmic decay is a valid approximation increases, cf. Fig. [3] 
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On a fixed point on that surface the decay is concave for q < q 
and convex for q > q, cf. Fig.0] For fixed q, the change from 
concave to convex is achieved by crossing the mentioned sur- 
face from above in the sense exemplified in Fig. [5] 

The coupled quantities share the leading asymptotic behav- 
ior of the density correlators. As a consequence of the factor- 
ization theorem of MCT, only the glass-form factors and the 
critical amplitudes h q are different for the coupled quantities 
ll33ll . The leading corrections imply a violation of a gener- 
alized factorization theorem. These are proportional to the 
correction amplitude K q . Since for large wave vectors, say 
q > 10, the quantities f q , h q \ and K s q are close to the ones for 
the coherent correlator, the approximation for the tagged par- 
ticle correlation functions § q {t) for these large q is the same 
as for <L(?). So the discussion for §i(t) is already exhausted 
by Fig.Q Not much could be gained from repeating the dis- 
cussion of the previous section for § q (t). 



IV. MEAN SQUARED DISPLACEMENT NEAR AN 
A4-SINGULARITY 




FIG. 6: Subdiffusive power law in the mean-squared displacement 
(MSD). The solutions for states 1, 2, and 3 in the inset are shown as 
full lines in the full panel together with the leading (dotted) and next- 
to-leading (dashed) approximation by Eq. 1 1 81 . The long horizontal 
line represents 6r* 2 = 0.01086, the short horizontal lines the correc- 
tions to the plateau, 6{r* 2 — r 2 ), cf. Eq. I19cl . The straight full lines 
show the power law (t/z) x , Eq. <21a> . with exponents x = 0.365, 
0.173 and 0.0878 for states n = 1, 2, 3. The filled symbols show 
the points where the solutions deviate by 5% from the leading-order 
power laws. The inset shows part of the glass-transition diagram for 
6 = 6* and a chain line where b2 = 0, cf. Eq. <20b> . (see text). 

According to Eq. j20al i. 5r 2 (?) is expected to exhibit power- 
law behavior around the plateau 6r£ 2 provided the term b2 
vanishes. The power-law exponent x is determined explicitly 
in Eq. J2 lbl by the localization length and the critical ampli- 
tude, which are 

r* = 0.04255 , ^ SD = 0.00405 1 . (24) 

The inset of Fig. [6] shows the line where bi from Eq. 



J20a> vanishes. This line is almost identical to the one for 
Z?2(24.2) = shown in Fig.[2]for the correlators § q (t). The 
MSD for three states on that line is shown in the full panel. It 
is described well by the approximation in Eq. Jl 8I >. For states 
n = 1, 2, 3, one, three and six decades are covered with devi- 
ations less than 5%, so the approximation yields a description 
of similar accuracy as for the correlation functions in Fig. [3] 
The leading result from Eq. describes the relaxation pro- 
portional to lnf (dotted) which always has negative curvature 
in the double-logarithmic representation and does not provide 
a valid description for n = 1 and 2. The reason for the qual- 
itative difference between the solution for the MSD and the 
leading logarithmic law is that the corrections proportional to 
^msd = — 1 -708 are large, Therefore, KmsdB 2 + B2 is never 
close to zero in the liquid regime except very close to the A4- 
singularity. This is seen for n = 3 in Fig.[6]where lnf develops 
a straightened decay around the plateau. 

The power law J2 It provides a different formulation of a 
leading order approximation and is shown in Fig.[6]as straight 
line for n = 1, 2, 3. For n = 1 this describes the MSD for more 
than a decade as indicated by the squares. For n = 2 three 
decades are covered and six decades of power-law behavior 
are identified for curve n = 3. So the accuracy is similar to the 
one provided by the approximation in next-to-leading order 
by Eq. Jl 81 . Both asymptotic descriptions fall on top of each 
other around the plateau and therefore corroborate that the re- 
formulation (j^Ua) is justified. The interpretation of the behav- 
ior of the MSD is then much simpler when considering the 
power laws instead of the logarithms of time. The decreasing 
slope of the relaxation when approaching the A4-singularity 
as in Fig. [5] is just the exponent x from Eq. J21bl which de- 
creases as B with the square-root of the separation parameter 
Ei, cf. Eq. Q. The same parameter B is the prefactor of the 
leading-order logarithmic decay in Eq. @. In that sense Fig. 
|5]is the analog of Fig. 

The term bi in Eq. J20at varies regularly in the separa- 
tion parameters £1 and £2, and bj. is positive above the line 
Z?2 = and negative below. Therefore, similar to the case for 
the correlators in the linear-log representation, in the double- 
logarithmic representation, the behavior of the MSD can be 
changed from convex to concave when crossing the line of 
vanishing bi- This is demonstrated for three states in Fig. 
State b is identical to the state n = 2 in Fig. [6] and obeys 
Z?2 = 0. The power law (f/x) x is shown as straight full line. 
The time scale x is matched at the plateau 6r* 2 . Moving to 
state c below the chain line, (r,cp) = (3.42,0.525), a relax- 
ation is obtained which clearly exhibits negative curvature and 
is consistent with the calculated value bi = —0.00735. The 
leading-order power law with exponent x = 0.285 fulfills a 
5%-deviation criterion for two decades which accidentally ex- 
tends to short times as the approximation crosses the solution 
twice. Reducing the allowed deviation to 4% would reduce 
that interval to less than a decade. If we include the term pro- 
portional to £>2 from Eq. J22bl and renormalize the exponent 



to x , Eq. ( I22at . the approximation agrees with the solution for 
three decades. It is obvious from a comparison with curve 1 
in Fig. [6] that the leading-order power law describes that solu- 
tion better than it describes the solution at state c in Fig.Qfor 
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FIG. 7: Concave and convex deviations from the power law, Eq. <2H 
in the MSD. Solutions for the states a, b, and c are shown as full lines, 
the approximation 12 1 1 as straight full lines for exponents x = 0.147, 
0.173, and 0.285, respectively. Filled symbols denote the 5% devia- 
tion of the solutions from the leading-order power law. For state b, 
the dashed line exhibits the corrected power law with x 1 = 0.155, Eq. 
<22a> . and the open triangle the 5% deviations of the solution from it. 
Dashed lines show the approximation by Eq. {22b) for a and c with 
b 2 = 0.00363 and -0.00735, respectively, and x 1 = 0.143 and 0.214. 
The open symbols mark the 5% deviations. The inset replots the one 
from Fig.|6|and shows by the crosses the state points a, b, and c. 



comparable values for X and the plateau correction . The de- 
viation to convex behavior is demonstrated by the dashed line 
at curve a, (r,(p) = (4.57,0.523). Again the range of validity 
is extended to earlier times but for later times no improvement 
can be found. 

In Fig. |6] the dashed line, which describes the next-to- 
leading order approximation of Eq. deviates from the 
leading order power law ( 12 1 at below the plateau where the 
range of validity for the power law extends to much smaller 
times than justified by its derivation. We also recognize that 
the exponent x overestimates the slope of the relaxation in 
Figs-EJandQ In Eq. d2 lbt only the term B from the leading 
order approximation is present. Taking into account the renor- 
malization of this prefactor to B — B\ in Eq. ( I22a> changes the 
exponent for state b from x — 0. 173 to x! = 0. 155. By compar- 
ing the full line for the leading result and the dashed line for 
the corrected one in Fig.Q we find that the range of applica- 
bility is shifted to later times by one decade and extended by 
two decades. The corrected power law is valid from t = 10 2 
to t — 5 • 10 6 and comparison to Fig.|6]shows that approxima- 
tion (|8} covers a similar range. The accidental extension to 
shorter times is removed. The approximation now covers the 
range also a naive power-law fit would yield. 

In summary, the correction amplitude Xmsd f° r the MSD 
does not vanish within the liquid regime. Therefore a log- 
arithmic relaxation law can be detected only for states very 
close to the singularity. However, there is a line of vanishing 
corrections for the logarithm of the MSD. Here a logarithmic 
relaxation can be observed and this describes a subdiffusive 
power law of the MSD. We can interpret Fig.0as the analog 



of Fig. |3 Some quadratic correction to a leading-order linear 
behavior can be set to zero on a surface in control-parameter 
space. Departing from that surface in opposite directions in- 
troduces either positive or negative corrections and the linear 
behavior is changed to convex or concave. 

V. A 3 -SINGULARITIES 




FIG. 8: Glass-form factors /° and amplitudes h° and h s °, and 
correction amplitudes K° and K q ° at the A3-singularity for 8 = 0.03. 
Line styles are the same as in Fig.Q The values at the Azj-singularity, 
fg (dotted), h* t (dash-dotted), and Kg (dash-dotted), are shown for 
comparison. The values for q = 24.2 and 45.0 are marked by dia- 
monds. The inset shows K„ for 4 < q < 11 for 8 = 8* (chain line), 
0.03 (full line) and 0.02 (dotted line). 

An A3 -singularity is not located on a liquid-glass-transition 
line but is the endpoint of a glass-glass-transition line yfl. 
The parameter /J3 is no longer vanishing and for 8 = 0.03 
we get /J3 = 0.109 and £ = 0.157. For this A3 -singularity the 
g-dependent amplitudes are shown in Fig. [8] No qualitative 
changes are obvious compared to the results shown in Fig.[T] 
for the A4. The smaller length scale 8 = 0.03 for the attractive 
well introduces a smaller localization length, and this implies 
the broader distributions in wave-vector space. So the trend 
seen when changing from the HSS to the A4-singularity of 
the SWS is continued when approaching A3-singularities at 
smaller 8. There are only two notable exceptions at smaller 
q. First, the value for K q at the position of the structure factor 
peak is minimal for the A 4 , - 1 .81 = K* < K° = - 1 .72. The 
inset shows this region enlarged for 8 = 8*, 0.03 and 0.02, 
demonstrating that Kg at the peak is again larger for the A3- 
singularity with smaller well width 0.02, where K q = —1.69. 
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Second, the zero-wave-vector limit of K s q is also smallest at the 
A4-singularity. The respective values for 8 = 8*, 0.03 and 0.02 
are —1.71, —1.64, and —1.62. Therefore, one experiences the 
strongest g-dependent corrections at the A4-singularity. 
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FIG. 9: Curves of vanishing quadratic correction in Eq. {8} for the 
A3 -singularity (O) at 8 = 0.03. The 8 = 0.03 cut through the glass- 
transition diagram is displayed by full lines. The various lines are 
shown in the same style as in Fig. [2| and labeled accordingly. The 
line Z>2 = 0, cf. Eq. I20bi . indicates the analogous line for the MSD, 
cf. inset of Fig.|6| 

Figure |9] shows the analog of Fig. [2] for a cut through the 
glass-transition diagram at 8 = 0.03. The lines 81 = and 
£2=0 for the A3 -singularity are obtained from a smooth trans- 
formation of the corresponding lines at the A4-singularity, and 
they appear in similar locations in the diagram. The line £2=0 
is again very close to the almost horizontal line of transitions. 
Just below, we find again the line where Bi{q) = when 
K CI = 0. However, this now represents q ps 57.5, cf. Fig. |8j 
which is a value almost twice as large as for the correspond- 
ing line in Fig. [2] For the wave vector q = 24.2 we find the 
line, where Z?2(24.2) = 0, completely in the glass state. Tak- 
ing the same value for the correction amplitude as for q = 24.2 
at the A 4 , K* ps -0.6, we obtain q = 45.0, cf. Fig. |S] and 
the line labeled accordingly in Fig. [9] Since the latter line 
comes close to the liquid-glass-transition line we take that as 
a reference and estimate the range of wave-vectors where the 
quadratic corrections can be put to zero in the liquid regime 
to 45 < q < 70. The lines where B2(q) — can be rather 
sensitive to ^-variation. This is demonstrated by the curve 
Z?2(46.2) = 0. Although the change in the wave vector is rela- 
tively small in comparison to q = 45.0, the values for K q differ 
by more than 20% for fixed q and induce a rotation of the line 
Bi(q) = by quite a significant angle. 

Having in mind the drastic changes in the lines where 
Bi{q) — 0, it may come with some surprise that the line for 
the MSD, where bi = 0, stays rather robust and accessible in 
the liquid regime as seen in Fig. |9] The variation in q for the 
amplitudes is reflected in changes of the localization lengths. 



For the A3-singularity at 8 = 0.03 we get 

r° = 0.0243 , h° MSD = 0.00136 . (25) 

From Eq. (I24> one gets r* jr° = 1 .75 and the square of the 
latter ratio, r* 2 /r° 2 ps 3, is the same as ^msd/' 2 msd- Since 
only the fraction liMSo/r 2 could introduce larger modifica- 
tions in Eq. j20aL the changes in b^ cancel approximately and 
the line specified by bi = experiences only minor deforma- 
tions when 8 is varied. The wave vector for which the lines 
Bi{q) = and b2 — Q are closest to each other, is q = 45.8 at 
the A3 -singularity for 8 = 0.03. 
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FIG. 10: Parameters for the asymptotic description at the A3- 
singularities of the SWS for varying 8. Panel A displays /J3 (A), 
Eq. O, and £ (0), Eq. {T^. The dashed curve shows the asymp- 
totic V8* — 8-law for the /J3. The localization length r° is shown in 
panel B. The correction amplitudes ^ M g D and the ratios ^msdAmSD 
are shown in panels C and D. 

To corroborate the finding for the MSD from the preced- 
ing paragraph, the parameters for the asymptotic description 
of the MSD at the A3 -singularities are shown in Fig.^| The 
/J3 vanish when we approach the A4-singularity. The decrease 
close to 8* is described asymptotically by a s quar e-root varia- 
tion, /J3 °= %/8* — 8, shown by the dashed line [33]. The small- 
ness of /J3 indicates that all the A3 -singularities are already in- 
fluenced by the proximity of the close-by A4-singularity. One 
can take advantage of this finding and conclude that the terms 
proportional to /J3 in Eq. (|8} are small. Moreover, one may ne- 
glect B3 and B4 in Eq. il lcl entirely without introducing large 
additional errors. The leading correction to the logarithmic 
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decay laws is then only quadratic also for the A3-singularities. 
Parameter X, varies regularly around a finite value at 8* but 
shares the variation of /J3 at 8* due to Eq. Jl 3I >. Panel B shows 
the decrease of the localization length at the A3 -singularity 
when 8 is reduced. A change of 40% in r° from 8 = 8* to 
8 = 0.03, cf. Eqs. J24l> and i25\ . is reflected in the broad- 
ening of the distributions in q seen in Figs. ^3 and |8] This 
broadening is responsible for the large variation in q when 
comparing Fig. [2] with Fig. [9] It was noted in the discussion 
of the inset of Fig.|S]that K q introduces the strongest correc- 
tions for the correlation functions at the A4-singularity. This 
is also true for the MSD as seen in panel C for /Tmsd which 
is largest in absolute value at the A4-singularity. The varia- 
tion in /^msd with 8 is however small and does not introduce 
significant changes to 02 in Eq. J20bl . The amplitude /imsd 
is the remaining parameter entering Eq. d20bi that could alter 
the location of the line b 2 = in the glass-transition diagram. 
We noted above that only the ratio ^msd/^ 2 nee ds to be con- 
sidered which is shown in panel D. From there one infers that 
the ratio varies only by less than 5%. We can conclude that 
the line of power law variation for the MSD stays in the liquid 
regime even when 8 is changed significantly. 
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FIG. 1 1 : Cut through the parameter space for the hard-core Yukawa 
system for 8 = 8*. Lines styles are the same as in Fig.[2| The wave 
vectors q = 15.0, 21.4, 27.0, and 28.2 are approximately equivalent 
to q = 7.0, 20.2, 24.2, and 27.0 in Fig.[5| respectively, after rescaling 
r by a factor of 2.98 (see text). 



VI. HARD CORE YUKAWA SYSTEM 

The A; -singularities occurring in MCT are topologically 
stable, smooth changes in the control parameters do not chal- 
lenge their existence. Therefore, the results for the SWS can 
be applied also to other potentials with a short-ranged attrac- 
tion. Nevertheless, the deformation of the potential might in- 
troduce changes large enough to be relevant for the detection 
of the higher-order singularities. Among several potentials the 
hard core Yukawa system (HCY) was found to differ by up to 
20% in certain properties at the A4-singularity from the SWS 
fl3ll . Since other potentials differ less we use that system as a 
second example for an A4-singularity. 

Figure ITTI shows the analog of Fig. [2] for the HCY. For a 
comparison, the A4-singularity in the SWS was mapped on top 
of the A4-singularity in the HCY by scaling in T with a factor 
of 2.98 and by a shift in 9 of 0.0065. The same transformation 
was applied to the lines where 62(9) = in the SWS. Fig. II II 
displays the B%(q) = lines for the HCY that come closest 
to the ones shown in Fig. Rafter the mapping. The correction 
amplitude K q for the HCY vanishes at q w 34, and the range in 
wave vector for which B 2 (q) = is lying in the liquid regime 
is shifted to higher wave vectors, 21 < q < 36 or —0.9 <K q < 
0.2, in comparison to the SWS. For q — 27.0 we get the line 
B 2 (q) = for the HCY that is closest to the line b 2 = for the 
MSD as compared to B 2 (24.2) = in the SWS. 



VII. CONCLUSION 

Logarithmic decay or, equivalently, 1// noise in the fluc- 
tuation spectra, can arise in a number of situations and is ex- 
plained by various approaches |36]. In the log-linear repre- 
sentation appropriate for the correlation functions, this decay 



exhibits a straight line. To discriminate the logarithmic decay 
laws originating from higher-order glass-transition singulari- 
ties within MCT 1 24] from other possible scenarios one needs 
criteria to distinguish one from the other. The theory makes 
specific predictions where in the control-parameter space the 
logarithmic decay is expected and how the corrections intro- 
duce deviations from that behavior. In this paper, the scenarios 
are discussed in quantitative detail for an example relevant for 
studies of colloidal dynamics, the square-well system (SWS). 
To proceed, specific cuts through the three-dimensional pa- 
rameter space are considered. Here, lines are identified where 
the corrections quadratic in the logarithm of time vanish for a 
chosen wave vector q, cf. Fig. |2] These lines emanate from 
the higher-order singularity and rotate clockwise around the 
higher-order singularity with increasing q. The correlation 
functions for states on these lines exhibit decays that are lin- 
ear in the logarithm of time for several orders of magnitude 
in time, cf. Fig. [3] In leading order, the slope of the decay is 
given by the square-root of the distance from the higher-order 
singularity, Eq. Q. The mean-squared displacement MSD 
displays a power law, Eq. ( 121 al l, that is valid on a similar line 
in the control-parameter space, cf. Fig. [6] The exponent x of 
this subdiffusive behavior is also decreasing with the square- 
root of the distance. Both the logarithmic decay and the power 
law are accessible in the liquid regime. The logarithmic decay 
is predicted for wave vectors q, which are equivalent to values 
of about three to four times the first peak of the static structure 
factor. 

In a semi-logarithmic representation for the correlation 
functions and a double-logarithmic plot of the MSD, char- 
acteristic convex and concave relaxation patterns are found 
when states are chosen that are off the specified lines, cf. 
Figs. [5] and Due to the variation of the correction ampli- 
tude K c/ in Fig.^ a similar variation from convex to concave 
behavior is introduced by changes in the wave vector at a fixed 
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point in control-parameter space, cf. Fig.|4] These deviations 
from logarithmic behavior provide a test for the clear identifi- 
cation of dynamical scenarios that are consistent with Eq. (|8) 
and hence originate from higher-order singularities. 

When the localization at the higher-order singularity is 
changed by either deforming the shape of the potential or by 
moving to A3-singularities at smaller ranges of the attraction, 
the logarithmic decay of the correlation functions is shifted 
to higher wave vectors. Whereas the difference between the 
SWS and the hard-core Yukawa system at the A4-singularity is 
modest, cf. Fig.|2]and Fig. ^2 the li nes of vanishing quadratic 
correction change drastically at the A3-singularity, cf. Fig.|9] 
In contrast, the line where the subdiffusive power law for the 
MSD is valid, is robust against changes of the well width and 
the potential shape, cf. Figs.|6][9] andll II 
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FIG. 12: Correlators for states a (panel a) and b (panel b) from Fig. 
Hfor wave vectors q = 4.2, 20.2, 24.2, 27.0, 32.2, and 36.2 from 
top to bottom. Full lines show the solutions of the MCT equations 
for the SWS, dashed lines the approximation by Eq. |8j. Triangles 
mark the 5% deviation of the correlator from the approximation for 
q = 27.0 and 24.2, respectively. The dotted vertical lines indicate the 
time scales x, the short horizontal lines the corrected plateau value 
fq + fq for q = 27.0 and 24.2, respectively. 

For comparing the solutions of the equations of motion, 
Eqs. Q, Q, and with the asymptotic expansions, Eqs. 
•ID, (1151 , and dl 81 . all parameters are calculated explicitly ex- 
cept the time scale X which is matched at the plateau. In an 
experiment or a computer simulation only the correlators are 
available directly. We show these in Fig. [21 for two states 
specified in Fig. [2] for different wave vectors. Since state b 
is closer to the Azt-singularity, the range of validity for the 
asymptotic approximation is larger than for state a. Espe- 
cially the extension of the linear- log decay at some specific 
wave vector increases when moving closer to the singularity. 
As noted in connection with Fig. |4] the range of validity for 
the approximation by Eq. ^ may vary with q. Partly for that 
reason a larger absolute curvature is attributed to the correla- 
tors by the approximation than a fit would do. A free fit could 
identify logarithmic behavior at state b for q = 20.2 from t ss 5 



to t 5 • 10 7 with a deviation of at most 5%. In addition, fit- 
ting the correlator for q = 24.2 also for t ^ 10 5 would yield 
positive curvature. Therefore, a free fit in that region of the 
control-parameter space tends to find the logarithmic decay 
at a somewhat lower wave vector than predicted by Eq. (jSJl. 
However, with a choice of the time scale x that is reasonably 
close to the theoretical value, the concave and convex decay 
patterns can still be identified unambiguously in the correla- 
tors without invoking additional assumptions. 
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FIG. 13: Comparison of f q and h q from the fit to the simulation of 
two different states 1 37] with the values ft and h* q for the SWS from 
Fig-0 For the comparison in the lower panel the theoretical values 
are multiplied by 0. 14. 

A recent molecular dynamics study of a binary mixture of 
square-well particles identifies a power law with x = 0.28 for 
the MSD over four decades and a related logarithmic decay 
of the correlation function at a wave vector q — 16.8 fyfa . A 
scenario similar to Fig.^l was found for the correlation func- 
tions: Upon increasing q, a change from concave to convex 
decay is observed. For a second state, faster decay with larger 
prefactors for the logarithmic decay is reported together with a 
larger exponent, x = 0.44, for the power law in the MSD. This 
finding is consistent with the assumption that this second state 
is further from the supposed higher-order singularity than the 
first state. Different from Fig. in the simulation 8 was 
changed to vary the distance while (p and T were kept fixed. 
The logarithmic decay was shifted to a higher wave vector for 
smaller 8 1 37] . This is consistent with the expectation that can 
be inferred from Figs.|5]and[5]by observing, e.g., the rotation 
of the line Z?2(24.2) = 0. The analysis of the simulation data 
allowed for a fit of the values for /* and h* 1 37] . These are 
shown in Fig. Altogether with the theoretical predictions for 
the SWS. The fitted parameters for both states almost fall on 
top of each other for f* . The amplitude h* is deduced from 
the simulation data only up to some overall factor. It can be 
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matched reasonably by a multiplication of the theoretical pre- 
diction for h* . The extension in q for the values obtained from 
MCT for the SWS are narrower, the width at half maximum 
for fq differs by 15%. A similar difference was observed for a 
binary mixture of hard spheres and agreement between theory 
and simulation could be improved by using the structure fac- 
tor from the simulation as input to the MCT calculations 1 38]. 
For the amplitude h* the locations of the maxima disagree by 
15% and the width is different by 25%. The deviations for 
q < 7 in both f q and h q can be attributed to the effects of mix- 
ing in. 

In summary, scenarios for logarithmic decay near higher- 
order glass-transition singularities are presented in this work. 



Some essential predictions are supported by the results of 
computer simulations. This should motivate further investi- 
gations in colloidal systems with short-ranged attraction. In 
particular the power-law behavior for the MSD including the 
deviations might be accessible to experiments. 
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